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ABSTRACT 

We present a demonstration of near real-time spacecraft astrometry with the VLBA. We detect 
the X-band downlink signal from Mars Reconnaissance Orbiter and Odyssey with the VLBA 
and transmit the data over the internet for correlation at the VLBA correlator in near real-time. 
Quasars near Mars in the plane of the sky are used as position references. In the demonstration 
we were able to obtain initial position measurements within about 15 minutes of the start of the 
observation. The measured positions differ from the projected ephemerides by a few milliarcsec- 
onds, and the repeatability of the measurement is better than 0.3 milliarcseconds as determined 
from measurements from multiple scans. We demonstrate that robust and repeatable offsets are 
obtained even when removing half of the antennas. These observations demonstrate the feasibility 
of astrometry with the VLBA with a low latency and sub-milliarcsecond repeatability. 

Subject headings: astronomical instrumentation — astronomical techniques 


1. Introduction 


The Very Long Baseline Array (VLBA) is an 
interferometer of ten 25 meter diameter anten¬ 
nas observing from 300 MHz to 90 GHz. Eight 
of the ten antennas are distributed throughout 
the continental United States, one is on Mauna 
Kea, H awaii and on e is on St. Croix, Virgin 
Islands (iNapier 1995ll . Long baselines (up to 
8611 km) make it possible to obtain an image 
angular resolution of 0.1 by 0.2 milliarcseconds 
at 8.4 GHz, and typical relative astrometric ac¬ 
curacies of about 100 to 200 microarcseconds 
for standard phase-referencing observ ations (e.g., 
Loinard et al.l 20071 : iMelis et~S^ 2014 ). These ca¬ 


pabilities allow the study of non-thermal astro¬ 
nomical phenomena at small angular scales and 
provid e precise astrometry, key to several applica - 


tions (iRomnev fc ReidI l2005t IZensus et al 


ppiica - 

Il998h . 


Precise astrometry has applications in spacecraft 
navigation as has been demonstrated by previous 
efforts by the National Radio Astronomy Obser¬ 


vatory (NRAO ) with the Spacecraft Navigation 


Pilot Project (iRomnev 2003; Lanvi et al. 20051: 


Martin-Mur et al.l 12006 ) and the Phoen ix Mars 
lander experiment ( Fomalont et al. 2010ll . These 


techniques have also been used to measure the 
mass of lapetus using data from a Cassini flyby 


(l.lacobson et al.N2006l ). observe Huygens as it en¬ 


tered Titan’s atmosphere to measure i ts position 
and help determine wind i n formation (iB ird et ^ 


2005 : Folkner et al. 200d Witasse et al. 2006ll . 


and to improve the ephe merides of Saturn b y us¬ 
ing Cassini observations ( Jones et al ] l20nll . All 
these projects used the standard VLBA observing 
mode in which data are recorded at the stations 
and shipped for correlation at the operations cen¬ 
ter in Socorro, NM. This process introduces a 
delay of about 2 weeks between observation and 
correlation of the data. 


The VLBA recently completed a sensitivity up¬ 
grade that included a new digital backend, a new 
Mark5C data recorder system and a new software 
correlator. The upgraded instrumentation at each 
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VLBA antenna includes a computer in the dig¬ 
ital data path which is primarily used for data 
routing to the recorder. This provides a conve¬ 
nient point to access the data allowing a reduc¬ 
tion in the latency of the data transfer from weeks 
to hours and thus enabling near real-time corre¬ 
lation. The correlation process has been simpli¬ 
fied and made mor e flexible by the D i stribu ted FX 
correlator (DiFX, IPeller et al. 2007 . 201lh which 
became the VLBA correlator in 2009 and has con¬ 
tinued to evolve in capability since then. These 
recent upgrades have opened the possibility of do¬ 
ing near real-time observations with the VLBA, 
which is a very attractive possibility for time criti¬ 
cal applications such as spacecraft navigation. To 
demonstrate the capabilities of the DiFX corre¬ 
lator for spacecraft navigation and exercise the 
use of the predicted ephemerides we performed 
observations of Voyager around t he time it was 
believed to enter interstellar space (iBurlaea et al 


2013t iKrimigis et al.ll2013t IStone et al.ll2013H . 


Here we present the results of a recent effort 
to demonstrate the feasibility of near real-time as¬ 
trometry of two spacecraft, Mars Reconnaissance 
Orbiter (MRO) and Odyssey, though passive ob¬ 
servations of their X-band downlink signals. We 
give a description of the observations including 
the preparatory work required to identify suitable 
phase calibrators for relative astrometry. We also 
describe the setup we use for near real-time data 
transfer over the internet, and the correlation pro¬ 
cess with DiFX using predicted ephemerides for 
each spacecraft. The data reduction process and 
results of these observations are described, includ¬ 
ing the latency, offsets with respect to predicted 
orbits, their accuracy and repeatability, and the 
redundancy and tolerance to missing antennas for 
the astrometry. We close the paper with a dis¬ 
cussion of some future developments required to 
make these observations standard and to improve 
the absolute astrometry. 


2. Observations 

This near real-time astrometry demonstration 
consists of observations at four epochs in 2013 
September for the MRO and Odyssey observa¬ 
tions and a preparatory observation of calibrators 
in 2013 August. 


2.1. Demonstration observations 

The preparatory observation, performed on 
2013 August 31, was a survey of possible phase 
calibrators near the position of Mars at the time 
of the near real-time astrometry demonstration. 
These data were used to determine the best and 
closest VLBA calibrator|i| within a set of candi¬ 
date sources less than 2 degrees from Mars at the 
time of the near real-time demonstrations. These 
calibrators should be bright enough to be detected 
in a short integration, and be nearly point like to 
minimize phase errors. Each one of the candidate 
sources was used to produce model images to use 
during the phase calibration of the demonstration 
observations. These observations were performed 
with standard latency with bandwidths of 1 MHz 
and 16 MHz. All the observations were reduced 
using standard methods for calibration and imag¬ 
ing in NRAO’s Astronomical and Image Process¬ 
ing System (AIPS jl. The 16 MHz bands were used 
to make the model images and measure integrated 
flux densities and angular sizes. From these re¬ 
sults we selected the most compact and brightest 
sources to be used as the calibrators, whose pa¬ 
rameters are given in Table [1] The 1 MHz images 
for the selected candidates were used to determine 
the minimum integration times required to obtain 
good phase solutions with the CALIB task on 
AlPS. In all cases we found good phase solutions 
even with 30 s of integration, which determines 
the minimum time required for phase calibrator 
integrations. 

The demonstration observations were per¬ 
formed on 2013 September 20, 21, 24 and 25. 
At each epoch we observed for 2 hours tracking 
Mars with both spacecraft in the same VLBA 
beam. Data were recorded in two 1 MHz bands in 
right circular polarization centered at the down¬ 
link frequency for each spacecraft, 8439.444446 
MHz for MRO and 8406.851853 MHz for Odyssey. 
To minimize data transmission, we use the lowest 
recording bandwidth available of 1 MHz. This pro¬ 
duces only 8 Mbps of data per antenna. We used 
simple phase-referenced observations by switching 
between calibrator (90 s of integration) and Mars 
(30 s of integration), with occasional observations 
of an astrometric check source. These astrometric 

^ http://www.vlba.nrao.edu/astro/calib/ 

^http://www.alps.nrao.edu 
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check sources have a known position and are used 
to assess the astrometric repeatability and accu¬ 
racy. We also observe a fringe finder source at the 
beginning, middle and end of the observations. 
For each one of the epochs the distance between 
source and calibrator was about 1.5, 1.4, 1.8 and 
1.2 degrees. A fixed Doppler correction was used 
which caused the central frequency to drift by 
about 100 kHz, mainly as a result of the space¬ 
craft orbital motion around Mars. This small drift 
is well within the observed 1 MHz band. 

A short video made by the NRAO outreach of¬ 
fice shows what happens in one of these demon¬ 
strations, and provides complementary informa¬ 
tion and a visual context to what is explained 
abov^l. 

Data transfer Most VLBA sites have limited 
network connectivity which is shared with mon¬ 
itor and control traffic from Socorro, NM. Two 
sites. Pie Town and Mauna Kea, can transfer at 
about 300 Mbps0 using 1 Gbps links funded by 
the US Naval Observatory (USNO) for their daily 
UTl—UTC observations. This opens the possibil¬ 
ity of wide-band measurements for a single base¬ 
line. All other sites can transmit about 1.1 Mbps 
on 1.4 Mbps links which are the limiting factor on 
the latency. Thus we can transmit about 3 GB of 
data per antenna in 6 hours, which corresponds to 
50 minutes of data at 8 Mbps and 12 seconds of 
data at the maximum recording rate of 2 Gbps. 

The near real-time data transfer is enabled by 
the XGube switct@, a Linux-PC based system sit¬ 
uated in the data path from the digital backends 
and the Mark5C recordei@. For data rates up to 
about 500 Mbps the data can be recorded onto a 
standard random-access filesystem on an internal 
hard drive on the computer. This is in contrast 
to the standard data path employing the Mark5C 
recorder where reading cannot be performed in 
parallel with writing. Files recorded to the inter¬ 
nal hard drive are transmitted over the internet us- 


^ http: //vimeo. com/80099626 

^Mbps: megabit per second, Gbps: gigabit per second 
®The XCube switch hardware is a product of the XCube 
Corporation http: //www. x3- c. com/ 

®Haystack Observatory, Mark 5 memo #57 
http://www.haystack.mit.edu/tech/vlbi/mark5/memo.html 


ing the Tsunami-UDP file transfer prograrrQ. The 
ability to continue transfers after interruption and 
its high performance over long-haul switched net¬ 
works leads to this program being commonly used 
for network transfer of Very Long Baseline Inter¬ 
ferometry data. 


Correlation There are several differences in the 
correlation of these observations compared to the 
vast majority of scientific observations. To accom¬ 
modate the spacecraft which are swiftly moving, 
the VLBA DiFX correlator is fed with the most 
recent ephemerides from the NASA’s Navigation 
and Ancillary Information Facility (NAIFjf|. The 
final positions determined do not depend in detail 
on the ephemeris; however, to prevent decorrela¬ 
tion (and hence loss of sensitivity) the ephemerides 
must predict the position sufficiently accurately 
(to a few arcseconds) as well as the rate of motion 
(to a few milliarcseconds per day). Because corre¬ 
lation occurred before the USNO published final 
values for Earth orientation parameters (EOF), 
predicted values had to be used. Since these pa¬ 
rameters are used to connect the interferometer to 
the celestial reference frame, errors in these pa¬ 
rameters result in errors in relative astrometry. 
Prediction errors for the EOP for a 5 day pre- 
diction horizo n have a maximum error of 5 mas 
I Malkin 2000ll . This results in a maximum offset 
error of about 87/ras for a I degree distance be¬ 
tween source and phase calibrator. The rms error 
of 2 mas for the same prediction horizon produces 
an rms offset error of about 35/xas. The VLBA 
DiFX correlator uses GSFG Gale version 9.1 as 
the VLBI delay model with additional corrections 
for wavefront curvature and the solar system grav¬ 
itational field to support near-field objects. The 
approximations made in the delay model start to 
break down at ~ 10® km distances. Two separate 
correlation passes were performed, one for each of 
the two spacecraft yielding two separate sets of 
corre lator output t hat are stored in FITS-IDI for¬ 
mat ( GreisenI 2011 ). 


2.2. Data calibration 

Calibration using AIPS starts with amplitude 
calibration including digital sampler bias correc- 


^http://tsunami-udp.sourceforge.net 
®NAIF: http://naif.jpl.nasa.gov 
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tion with the ACCOR task. During the near real¬ 
time observations no system temperature data are 
available until the end of the observing run, so 
system temperature calibration is skipped for the 
first stage of the data analysis. The data pre¬ 
sented here include the system temperature cali¬ 
bration. This does not affect astrometry as it sim¬ 
ply changes the absolute scale of the visibilities. 
Polarization parallactic angle rotation corrections 
are performed using the CLCOR task. Instrumen¬ 
tal delay corrections are done by fringe fitting on a 
bright calibrator at the beginning of the observa¬ 
tion using the FRING task. Instead of a bandpass 
calibration we simply exclude the 10% outer chan¬ 
nels, out of a total of 200, that have the largest 
gain variation. Phase calibration uses the nearby 
calibrator sources we modeled in the preparatory 
observations. For this last step we use the CALIB 
task. The phase calibration step registers the po¬ 
sition of the spacecraft to the ICRF. Examples of 
calibrated bandpasses for the MRO and Odyssey 
are shown in Figure [1] 

2.3. Imaging and astrometry 

Imaging is done using the IMAGR task in 
AIPS. An image is produced for each 30 s scan first 
using all the available antennas but later discard¬ 
ing any problematic ones to improve the image 
quality. The images are 1024x1024 pixels, with 
0.2 mas pixels. For MRO we use all the chan¬ 
nels, but for Odyssey only 60 channels centered 
on the peak are used. These channels are not the 
same for all scans due to Doppler shift that we are 
not compensating. Example images for a single si¬ 
multaneous observation of MRO and Odyssey on 
a single TS030E scan are shown in Figure [2j In 
the figure it can be seen that the artifacts in both 
images are almost identical due to correlated cal¬ 
ibration errors arising from the close proximity of 
the sources, which are separated by at most 5 arc- 
seconds over the different epochs and scans. This 
implies that our ability to measure the relative 
spacecraft positions is much better than their ab¬ 
solute ones (see the discussion in Section [3]). 

Astrometry is performed directly in the images 
by fitting a Gaussian model component with the 
JMFIT task which provides the offsets with re¬ 
spect to the phase center used for correlation. 


3. Results 

The main goal of these observations was to 
demonstrate the feasibility of near real-time as¬ 
trometry and not to provide optimal results. With 
this limitation in mind it is still useful to look 
at the characteristics of the results that this low- 
bandwidth observations calibrated with a simple 
scheme can provide. 

3.1. Repeatability and accuracy 

We estimate the repeatability of the astrometry 
by obtaining the offsets between the center of the 
fitted Gaussian component and the phase center 
for each 30 s scan on each spacecraft. The offsets 
are Aa in R.A. and A6 in deck Since each offset 
is measured using non-overlapping data they pro¬ 
vide independent measurements of the spacecraft 
position. The error budget for each measured off¬ 
set will contain contributions from the ephemeris, 
calibrator position errors and calibration transfer 
errors. A summary of the repeatability results for 
the MRO is shown in Tableland for the Odyssey 
in Table [3] 

From these results we can see that the measured 
spacecraft position was offset from the ephemeris 
by different amounts on each observation, with 
values up to a couple of milliarcseconds. The scat¬ 
ter in the measured offsets is less than about 0.3 
mas and in many cases less that 0.1 mas. 

We estimate the accuracy of the astrometry by 
using the astrometric check source observations, 
which are observed in the same way as the space¬ 
craft, but for a source of well determined position. 
A summary of the offset between the catalog and 
observed position of the calibrators is shown in Ta¬ 
ble m This offset can be as high as about 3 mas. 
This result indicates that the limiting factor for 
the accuracy is the quasar grid. This is a problem 
that could be minimized by obtaining astrometric 
observations on the calibrator sources. 

We also estimate the repeatability of the rel¬ 
ative offset between the spacecraft by taking the 
difference between estimated offsets for simultane¬ 
ous scans in Odyssey and MRO. As pointed out in 
Figure [2] and associated discussion, these relative 
offset should have a lower scatter than each indi¬ 
vidual spacecraft offset. The standard deviation 
in the offset difference is less than 0.055 mas in all 
cases, and is much better in most cases. See Table 
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[5] for a summary. 

3.2. Effect of reduced number of antennas 

The ten VLBA antenna^ can provide a highly 
robust system due to redundancy, which is key for 
applications requiring near real-time astrometry. 
In order to quantify the robustness of our results, 
we tested the effects in the astrometry of the loss 
of one or several antennas. This is tested by ex¬ 
cluding sets of antennas from a single TS030C scan 
and repeating the same astrometric measurements 
performed on the complete data for each scan. 
Only 9 antennas were used on the complete data 
set because KP had sensitivity problems at the 
time of observations. Omitting a single antenna 
(i.e., using 8 antennas) resulted in good quality 
images for the 9 cases we tested with rms offsets 
with respect to the original image of 35 /ras in a 
and 70 fias in <5. The maximum excursion was 
about 130 /ras, measured when removing SC. The 
second largest excursion, about 100 yitas, was ob¬ 
tained when removing BR. We tested by omitting 
selected pairs of antennas (i.e., using 7 antennas), 
and measured maximum excursions of about 200 
^as. Even larger offsets (up to 600 fias) are ob¬ 
served when including only 5 antennas (BR, FD, 
OV, LA, PT). These differences can be reduced by 
the use of improved calibration techniques and by 
increasing the time on source. 

One should notice that LA experienced tech¬ 
nical difficulties (excluded from TS030D/E). For 
TS030D, MK experienced high system tempera¬ 
tures during the second half of the epoch and was 
excluded. In the same epoch SC and HN were 
observing at elevations lower than 20 degrees and 
were also excluded. This resulted in TS030D ef¬ 
fectively using 5 antennas, and producing lower 
quality results than the other epochs. 

4. Summary and discussion 

We presented the results of a demonstration of 
near real-time astrometry for the MRO and the 
Odyssey performed at four epochs during 2013 
September. The feasibility of this observing mode 
was successfully demonstrated, with first detec¬ 
tion of the fringe calibrator in less than 5 min- 


®We refer to them by their code names: BR, FD, HN, KP, 
LA, MK, NL, OV, PT and SC. 


utes after the start of the observations and first 
spacecraft offsets within 15 minutes of the start of 
the observations. The results are promising, even 
with the simple observation, calibration, imaging 
and astrometric procedures used, indicating that 
significant improvements can be made by refining 
these procedures. The measured positions differ 
from the projected ephemerides by a few milliarc- 
seconds, and the repeatability of the measurement 
within a single epoch is better than 0.3 milliarc- 
seconds as determined from measurements from 
multiple scans. 

The technique is promising for spacecraft above 
—30 degrees of declination, beyond which calibra¬ 
tion and mutual visibility of the constituent anten¬ 
nas become increasingly difficult to handle. The 
limited network connectivity limits the choice of 
calibrators to bright sources of at least 50 mJy for 
the current 1 MHz bandwidth. The image plane 
analysis we performed requires a good correlator 
model. Total delays are obtained as part of the 
data analysis and are readably obtainable. 

Overall the demonstration was successful, pro¬ 
viding results with good accuracy, repeatability 
and tolerance to missing antennas. Straightfor¬ 
ward improvements to this technique can make it 
a valuable addition for spacecraft navigation. 

W.M. thanks Amy Mioduszewski for useful 
conversations about phase-referenced observations 
with the VLBA. The National Radio Astronomy 
Observatory is a facility of the National Science 
Foundation operated under cooperative agreement 
by Associated Universities, Inc. This work made 
use of the Swinburne University of Technology 
software correlator, developed as part of the Aus¬ 
tralian Major National Research Facilities Pro¬ 
gramme and operated under licence. We thank 
the referee for comments that improved the pre¬ 
sentation of the paper. 
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Table 1 

Properties of the calibrator sources. 


Source 

RA 

hh:mm:ss 

Dec 

dd:mm:ss 

RA error 

mas 

Dec error 

mas 

Size®^ 

mas Xmas 

Peak flux 
mjy/beam 

Note*’ 

Reference‘s 

J0927+3902 

09:27:03.0139380 

39:02:20.851770 

0.130 

0.100 

0.74x0.42 

6896±19 

Fringe 

GSFC 

J0908+1609 

09:08:55.9253500 

16:09:54.763880 

0.553 

1.393 

0.36x0.00 

195±1 

Phase/Check 

GSFC 

J0909+1928 

09:09:37.4434660 

19:28:08.291220 

0.140 

0.190 

0.74x0.00 

195±1 

Check 

PETR 

J0925+1658 

09:25:49.9644620 

16:58:12.203380 

1.472 

2.367 

0.08x0.00 

148±1 

Phase 

GSFC 


^This is the deconvolved source size major and minor axes. A size of zero indicates the source is unresolved in this direction. 
^ “Fringe” is for fringe fitting calibrator. “Phase” is for phase calibrator. “Check” is for astrometric check source. 

°GSFC: lMa et all {200^1 

PETR: L. Petrov, solution rfcJ2012b (unpublished, available on the Web at http://astrogeo.org/vlbi/solutions/rfc_2012b) 


Table 2 

Summary of the repeatability of the offset measurements for MRO. 


Obs. 

Date 

Scans 

[count] 

Aa 

[mas] 

[mas] 

^ Aa 

[mas] 

CTA5 

[mas] 

TS030B 

Sep 20 

8 

0.82 

0.12 

0.12 

0.15 

TS030C 

Sep 21 

12 

0.76 

0.06 

0.06 

0.13 

TS030D 

Sep 24 

9 

2.95 

-1.31 

0.23 

0.18 

TS030E 

Sep 25 

17 

2.42 

-1.47 

0.14 

0.06 


Table 3 

Summary of the repeatability of the offset measurements for Odyssey. 


Obs. 

Date 

Scans 

[count] 

Aa 

[mas] 

[mas] 

O' Aa 

[mas] 

(TAS 

[mas] 

TS030B 

Sep 20 

4 

0.38 

0.33 

0.10 

0.07 

TS030C 

Sep 21 

6 

0.26 

0.23 

0.19 

0.27 

TS030D 

Sep 24 

10 

2.22 

-1.08 

0.23 

0.15 

TS030E 

Sep 25 

12 

1.88 

-1.39 

0.09 

0.13 


Table 4 

Summary of the offset between catalog and observed calibrator position. 


Obs. 

Calibrator 

Check source 

Separation 

[deg] 

Aa 

[mas] 

A6 

[mas] 

TS030B 

J0908+1609 

J0909+1928 

3.3 

0.19 

-1.14 

TS030C 

J0908+1609 

J0909+1928 

3.3 

0.23 

-1.24 

TS030D 

J0925+1658 

J0908+1609 

4.1 

2.93 

-1.54 

TS030E 

J0925+1658 

J0908+1609 

4.1 

1.78 

-1.78 
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Table 5 

Summary of the repeatability of the offset measurements difference for Odyssey and 

MRO. 


Obs. 

Date 

Scans 

[count] 

Aa 

[mas] 

A,5 

[mas] 

[mas] 

CTAS 

[mas] 

TS030B 

Sep 20 

2 

-0.170 

0.163 

0.027 

0.016 

TS030C 

Sep 21 

4 

-0.651 

0.243 

0.054 

0.045 

TS030D 

Sep 24 

8 

-0.711 

0.219 

0.025 

0.009 

TS030E 

Sep 25 

9 

-0.405 

0.020 

0.007 

0.055 





Channels 

Lower frame: AmplJy Top frame: Phas deg 

Vector averaged cros^power spectrum Several baselines displayed 

Tlmerange: 0CV13:1627 to 0(V13:16:51 


Fig. 1.— Example of calibrated bandpass for MRO (right panels) and Odyssey (left panels). Each panel 
has phase in degrees on top and amplitude on the bottom, in kilo Jansky for MRO and Jansky for Odyssey, 
both as function of channel number. The upper panels are for the NL-SC baseline and the lower ones for 
the OV-PT baseline. 
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Plot file version 2 created 11-DEC-2013 20:44:16 
GREY:MRO IPOL 8439.747 MHZ MRO.ICL001.12 
CONT:MRO IPOL 8439.747 MHZ MRO.ICL001.12 

0 SO 100 150 200 



■10 I 


10 5 0 -5 -10 

MIIIIArc seconds 

Center at RA 09 22 09.5381360 DEC 16 41 20.031762 
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Fig. 2.— Example of calibrated image for MRO (left panel) and Odyssey (right panel), for a simultaneous 
observation of both spacecraft on a single TS030E scan. Coordinates represent the offset with respect to the 
phase center used for correlation. The MRO offsets are Act = 2.34 mas in right ascension and AS = —1.54 
mas in declination, with a signal to noise of 33.4. For the Odyssey the offsets are Act = 1.93 mas and 
AS = —1.57 mas, with a signal to noise of 32.8. The artifacts in both images are almost identical due to 
correlated calibration errors arising from the close proximity of the sources, which are separated by at most 
5 arcseconds over the different epochs and scans. 
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